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Context. Little is known about individual molecular clouds at high redshifts. The gravitational lens toward B0218+357 
offers the unique possibility to study cool moderately dense gas with high sensitivity and angular resolution in a cloud 
that existed half a Hubble time ago. 
■ Aims. This non-CO molecular multi-level study of a significantly redshifted cloud aims at a better definition of the 

physical properties of molecular gas in a kind of interstellar environment that is rarely studied in the Galaxy or in the 
nearby extragalactic space. 

Methods. Observations of the radio continuum and six formaldehyde (H2CO) lines were carried out with the VLA, the 
Plateau de Bure interferometer and the Effelsberg 100-m telescope. 

Results. Three radio continuum maps indicate a flux density ratio between the two main images, A and B, of ~3.4±0.2. 
Within the errors the ratio is the same at 8.6, 14.1, and 43 GHz. The loi — 0oo line of para-H2CO is shown to absorb the 
continuum of image A. Large Velocity Gradient radiative transfer calculations are performed to reproduce the optical 
depths of the observed two cm- wave "K-doublet" and four mm- wave rotational lines. These calculations also account for 
a likely frequency-dependent continuum cloud coverage. Confirming the diffuse nature of the cloud, an n(H.2) density of 
1— 1, <1000cm -3 is derived, with the best fit suggesting nffb) ~ 200 cm -3 . The H2CO column density of the main velocity 

component is ~5xl0 13 cm -2 , to which about 7.5X10 1 cm -2 has to be added to also account for a weaker feature on the 
blue side, 13kms _1 apart. iV(H2CO)/JV(NH3) ~ 0.6, which is four times less than the average ratio obtained from a 
small number of local diffuse (galactic) clouds seen in absorption. The ortho-to-para H2CO abundance ratio is 2.0-3.0, 
which is consistent with the kinetic temperature of the molecular gas associated with the lens of B0218+357. With the 
gas kinetic temperature and density known, it is found that optically thin transitions of CS, HCN, HNC, HCO + , and 
C**) ' N2H + (but not CO) will provide excellent probes of the cosmic microwave background at redshift 2=0.68. 
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'. 1. Introduction CS (Combes e t al. 11597)1 . OH (Kanekar et al. 120051) . NH 3 

^ ' n™,r, ■ r i • • i i (Henkel et al. 120051) and tentatively also LiH (Combes & 

B0218+357 is one of the most interesting gravitational lens wiklind ngggp , xhese absorp ti n lines originate from a 

systems studied to date. The background source showing diffuse molecular cloud of Wtic tem p era ture T kin ~55K 

a red featureless optical spectrum (Stickel & Kuhr [Mp (Henkd et &1 that ig located al the lme _ of _ sight 

is believed to be a BL Lac object The prominent radio to component A (Menten & Re i d [Tg0J|). The cloud param- 

source has a redshift of z~0.944±0.002 (Cohen et al. eters appear to be quite peculiar when compared with the 

and is lensed by an mtervenmg almost face-on spiral (for ties of clouds in the Galaxy or in nearby extraga i ac _ 

an image, see York et al. WHB mto two compact images, tk mg (Henkel ^ ai m ™ 

A and B, and afamt Emstem rmg (e.g., O Dea et al. [19921 gravitational lens of the B021 8+357 svstem nro 
Patnaik et al. [1993}. The flux density of source A is about . , ine gravitational iens ot the B0218+35 1 system pro- 
three times higher than that of source B; their separation vldes , a um( l ue laboratory for detailed molecular cloud re- 
as well as the angular diameter of the ring are with ~0'/3 f arch at intermediate redshftts. From the standard model 
exceptionally small. (A-cosmology with iJ =71kms Mpc \ fi m =0.27 and 
t> i. i M ftMi i 4. j i a i- j-- f f2A=0.73: Spergel et al. 2003) we obtain a light travel time 
Browne et al. (|1993p reported optical detections oi nar- 71 2 G r a luminosit distance of 4 1 G c a comovin 
row atomic absorption and emission lines from the lens at ,.' ,^ r ' a un } ln( i S l-,'^ 1S ance ' Pp' a comovm g 
redshift z=0.6847. H I absorption was observed by Carilli et ™ ^ distan f • ° 5 ^ PC ^ T ^ ^ 8126 -i 3 ^ 6 ° 
al. (HMD and Kanekar et al. J2JJ0J. The first molecular ab- L5 Gpc ' resultm S m a hnea r scale of 7.1 pc mas . 
sorption lines from t he int ervening system were detected by In the following we present three radio continuum maps 
Wiklind & Combes (flOPS]) . Identified molecu lar spe cies are 01 B0218+357 as well as the first non-CO molecular multi- 
CO, HCN, and HCO+ (Wiklind & Combes [TOOSD . H 2 CO l eve l stud y of a significantly redshifted cloud. The species of 
(Menten & Reid lT996jh H 2 (Combes & Wiklind fTMT)) . choice is formaldehyde (H 2 CO). With the detection of five 

H2CO transitions complemented by an additional H2CO 

Send offprint requests to: N. Jethava, e-mail: spectrum from Menten & Reid (|l!)!)(i|). we have a total of 

nj ethavaCsmpif r-bonn . mpg . de six lines (see also Darling & Wiklind I2005|) . These and the 
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continuum data can be used to further elucidate the highly 
peculiar physical cloud parameters (Henkel et al. 2005]), to 
constrain the continuum morphology and to estimate an 
ortho- to para-H2CO abundance ratio that has the poten- 
tial to provide clues to the history and chemistry of the 
cloud. 



2. Observations 

Table [T] and Fig.[T] summarize the observational parame- 
ters of the H 2 CO line measurements, including the 2 1:L — 2 12 
line reported by Menten & Reid ()1996|) . Observations 
were carried out with the Very Large Array (VLA) of 
the National Radio Astronomy Observatory (NRACQ) , the 
Plateau de Bure Interferometer (PdBI) of the Institute 
de Radioastronomic Millimetrique (IRAIvfl) and the 100- 
m Effelsberg telescope of the Max-Planck-Institut fur 
Radioastronomie (MPIfR0). 



2.1. VLA 

The VLA observations were carried out in the A- 
configuration at 8.6 and 14.1 GHz (X- and U-band; see 
Menten & Reid [TUMI) and at 43.2 GHz (Q-band; projects 
AM0494 and AM0545), employing 25, 25, and 13 25-m an- 
tennas, respectively. Measurements at the two lower fre- 
quencies, already presented by Menten & Reid 



(|1996|) . were 

re-reduced (note the slight difference in angular resolution 
of the X-band data between those given in our Table Q] and 
that of Menten & Reid HMrJ]) . In the X- and U-band, the 
total bandwidth of 6.25 MHz was split into 64 channels. 
In the Q-band, the bandwidth was 12.5 MHz, consisting 
of 32 channels. The phase center was located at aj2ooo = 
02 h 21 m 05!4733 and <5j 2 ooo = +35°56'13'.'791. Phase and 
bandpass calibration were achieved by observations toward 
the standard VLA calibrator 3C 84. Absolute flux calibra- 
tion was obtained by measuring 3C 286, for which flux den- 
sities of 5.11, 3.54 and 1.46 Jy were adopted at 8.6, 14.1 and 
43.2 GHz, respectively (Ott et al. [TgM|) . 

The data were calibrated and reduced according to the 
standard NRAO Astronomical Image Processing system 
(AIPS). Using the AIPS task UVLIN, we removed the con- 
tribution of the continuum in the uv-plane to obtain the 
spectral line data that were imaged and CLEANed. The 
B0218+357 continuum source is strong (see Sect. 3.1) so 
that several iterations of self calibration could be performed 
to improve the image quality and to determine the phase 
and amplitude corrections which were applied to each spec- 
tral line channel. The data were converted into CLASS for- 
mat to perform the final analysis. 

2.2. PdBI 

The PdBI observations were carried out in the 5C2- 
array configuration (project G048) at frequencies of 83.601, 



1 The NRAO is a facility of the National Science Foundation 
operated under cooperative agreement with Associated 
Universities, Inc. 

2 IRAM is supported by INSU/CNRS (France), the MPG 
(Germany) and the IGN (Spain). 

3 Based on observations with the 100-m telescope of 
the MPIfR (Max-Planck-Institut fur Radioastronomie) at 
Effelsberg, Germany. 
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Fig. 1. H 2 CO energy level diagram up to 40 K above the 
ground state. Note the "K-doublets" in the K a =l ladder, 
that are the result of the slight asymmetry of the molecule. 
For the six observed transitions that are marked by arrows 
(two J=2-l ortho-H2CO rotational transitions, one J=l-0 
and one J=2-l para-H2CO rotational transition and the 
J=l and 2 K-doublet ortho-H 2 CO transitions), rest fre- 
quencies are given. 



86.429 and 89.335 GHz (see Table □]). The array was 
equipped with five 15-m antennas. At 83 GHz, amplitude 
calibration was obtained using the standard PdBI calibra- 
tors 0415+379 and 0234+285 and adopting flux densities of 
6.28 and 1.36 Jy. At 86 GHz, 3C 454.4, 3C84 and 0415+379 
with fluxes of 6.14, 5.38, and 4.1Jy were used. Amplitude 
calibration at 89 GHz was obtained measuring 2230+114, 
3C84 and adopting fluxes of 7.97 and 4.64 Jy, respectively. 
For the phases, the standard PdBI calibrator 3C 84 was 
observed. The data were reduced using "Continuum and 
Line Interferometer Calibration" (CLIC) standard proce- 
dures. The calibration tables were fed into MAPPING to 
subtract the continuum levels and to convert the data into 
CLASS format. Since B0218+357 is a strong continuum 
source (Sects. 1 and 3.1), we could perform self-calibration, 
which significantly reduced formal phase and amplitude un- 
certainties. 



2.3. Effelsberg 

Observations with the 100-m telescope at Effelsberg were 
carried out in September 1997 using a dual channel 9 cm 
wavelength HEMT receiver. These measurements were 
done in a position switching mode with a two level auto- 
correlation spectrometer containing 2x512 channels with a 
bandwidth of 6.25 MHz. Pointing corrections could be de- 
termined toward B0218+357 itself and were found to be 
smaller than 10". Calibration was obtained from contin- 
uum cross-scans toward 3C286. 
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Molecule 


Line 


Frequency 


Telescope 


Date 


Channel 


FWHM 


Position 


rms 






Rest Redshifted 




of 


width 


beam 


angle 














observation 




size 










(GHz) 








(km s" 1 ) 


(arcsec) 


(degrees) 


(mjy) 


H 2 CO(0) 


1 10— In 


4.830 


2.867 


Effelsberg 


05-Sep-97 


1.6 


250 




1 


H 2 CO(0) 


2n— 2i 2 


14.488 


8.600 


VLA-A 


24-Aug-95 


1.7 


0.26x0.24 


27 


2 


H 2 CO(P) 


loi - Ooo 


72.838 


43.236 


VLA-A 


02-Dec-96 


3.1 


0.05x0.04 


16 


10 


H 2 CO(0) 


2i 2 — In 


140.840 


83.601 


PdBI-5C2 


09-Nov-96 


0.25 


6.03x5.29 


11 


36 


H 2 CO(P) 


2() 2 — loi 


145.603 


86.429 


PdBI-5C2 


23-Mar-97 


0.25 


5.83x5.06 


15 


31 


H 2 CO(0) 


2n — lio 


150.498 


89.335 


PdBI-5C2 


05-Aug-97 


0.25 


5.64x4.39 


14 


35 



a) Col. 1: H 2 CO species (O: ortho; P: para). Col. 2: Quantum numbers of the respective transition (JxaKc). Cols. 3 and 4: Rest and 
redshifted frequencies. Cols. 5-7: Observatory (VLA: Very Large Array; PdBI: Plateau de Bure Interferometer), epoch and channel 
width. Cols. 8 and 9: Full Width to Half Power (FWHP) beam size and major axis beam orientation. Col. 10: la noise level. 



Table 2. Continuum flux densities") 



Frequency 


Epoch 


Total flux 


A/B 






density 




(GHz) 




(Jy) 




2.9 


05-Sep-97 


0.69±0.06 




8.6 


24-Aug-95 


1.30±0.13 


3.3±0.13 


14.1 


24-Aug-95 


1.18±0.12 


3.3±0.14 


43.2 


02-Dec-96 


0.34±0.03 


3.5±0.22 


83.6 


09-Nov-96 


0.58±0.09 




86.4 


23-Mar-97 


0.48±0.07 




89.3 


05-Aug-97 


0.52±0.08 





a) Col. 1: Observed frequencies; to compare with the rest fre- 
quencies given in Tables [T] and [3] multiply by 1.68. For details 
of the 8.6 and 14.1 GHz masurements, see also Menten & Reid 
(1996). Col. 4: Flux density ratios between the two main lensed 
images of the background source. 

3. Results 

3.1. Continuum 

Menten & Reid (|1996|) presented VLA A-array continuum 
maps of B0218+357 at 8.6 and 14.1 GHz. We have re- 
reduced their data along with new A-array observations at 
43.2 GHz. Fig. [5] presents the three continuum images. The 
upper panel displays the 43.2 GHz map with the two main 
radio continuum components A and B being well separated. 
The two components and traces of the ring (e.g. Mittal et 
al. l2006[) are seen in the central panel. These main features 
are less well separated in the lower panel, where the an- 
gular resolution is not quite as good. Total flux densities 
including the Effelsberg and PdBI data as well as the VLA 
A-array flux density ratios between images A and B are 
given in Tabled 

3.2. Formaldehyde absorption toward image A 

All five formaldehyde lines we searched for are detected. 
Including the re-reduced 2n — 2i2 line reported by Menten 
& Reid (|1996j) . Fig. [3] shows all H 2 CO lines measured to 
date toward B0218+357. The upper two panels show "K- 
doublet" lines. Rotational lines are shown in the lower four 
panels. In order to minimize systematic calibration errors 
and to focus on physical parameters, line-to-continuum ra- 



tios are given that provide a direct estimate of optical 
depths. Such ratios should also be free of problems related 
to variability of the calibrator, to variability of B0218+357 
itself (e.g., Biggs et al. l2001a[l2001bj) or to elevation depen- 
dent gain variations of the respective telescope(s). 

Two of the measured H2CO transitions belong to the 
para and four to the ortho species (Fig. [3]). Systematic dif- 
ferences between para and ortho-H2CO lines are not ap- 
parent. Instead, there are obvious differences between the 
two K-doublet profiles measured at cm-wavelengths and 
the four rotational transitions observed at mm-wavelengths. 
The former show line-to-continuum ratios that are far be- 
low unity, while the latter are characterized by much higher 
line-to-continuum ratios. And while the former only show 
one velocity component, centered at about +6kms _1 , the 
latter show a second component at about -7kms _1 . It is 
not clear whether this difference is intrinsic to the source 
or whether it is merely an effect of limited signal-to-noise 
ratios in the K-doublet lines. Line parameters obtained by 
Gaussian fits are given in Table [3l 

Menten & Reid (fl"996|) found that absorpt ion in the 
2n — 2i2 line of ortho-H 2 CO is confined to source A. The 
loi — Ooo para-H 2 CO transition, that was taken with even 
higher angular resolution (see Table [1]), yields the same re- 
sult. There is no hint of line absorption toward source B 
with a 3.1 km s -1 channel 3er upper limit to the line-to- 
continuum ratio of 0.39 (see also Muller et al. I2007p . 

4. Discussion 

4.1. Radio continuum maps and spectral indices 

Table [2] indicates that the flux density of B0218+357 
reaches a peak near 10 GHz. However, B0218+357 is a vari- 
able source (e.g., Biggs et al. I2001b[) and our data were 
taken during a time interval of two years. In order to obtain 
a reliable radio continuum spectrum, quasi-simultaneous 
observations at a number of frequencies are mandatory. 
Mittal et al. (2006) presented such observations at five fre- 
quencies ranging from 1.65 to 15.35 GHz. The combined 
flux density of sources A and B (excluding the Einstein 
ring) reaches indeed a peak near 10 GHz (0.8-0.9 Jy; Mittal 
et al. 2007). Our total fluxes (Table \2§ include the ring and 
are thus higher. 

The frequency range covered here by the continuum 
maps is larger than that of Mittal et al. (2006) . Expecting 
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Table 3. H 2 CO line parameters ) 



Line 


Species 


Rest 


J'(T L /T c )dV 




AVi/2 


|T L |/T C 


Sa/ Stot 


Optical 








frequency 










depth 








(GHz) 




(kms^ 1 ) 










lio 


-In 





4.82 


-0.073±0.006 


5.3±0.2 


5.1±0.5 


0.014±0.002 


0.56 


0.025±0.006 


2n 


-2l2 





14.48 


-0.038±0.004 


5.1±0.4 


8.7±0.9 


0.004±0.0006 


0.70 


0.006±0.002 


loi 


-Ooo 


P 


72.83 


-1.154±0.097 


5.7±0.2 


5.3±0.5 


0.217±0.028 


0.77 


0.331±0.047 










-0.175±0.066 


-6.9±0.4 


2.7±0.5 


0.065±0.027 




0.088±0.036 


2l2 


-In 





140.83 


-1.721±0.063 


5.8±0.1 


5.7±0.3 


0.300±0.018 


0.77 


0.493±0.051 










-0.304±0.046 


-6.9±0.2 


3.1±0.6 


0.098±0.024 




0.136±0.035 


2o2 


-loi 


P 


145.60 


-1.047±0.052 


6.0±0.1 


5.1±0.4 


0.206±0.018 


0.77 


0.311±0.034 










-0.105±0.029 


-7.3±0.3 


1.9±0.7 


0.056±0.026 




0.076±0.034 


2n 


-lio 





150.49 


-2.458±0.102 


5.8±0.1 


5.4±0.3 


0.454±0.028 


0.77 


0.890±0.071 










-0.363±0.056 


-7.4±0.2 


2.2±0.4 


0.164±0.037 




0.240±0.053 



a) Cols. 1 and 2: H2CO transition (O: ortho; P: para). Col. 3: Rest frequency. Cols. 4-7: Integrated line-to-continuum ratios, 
heliocentric velocity with respect to 2=0.68466, linewidth, and peak line-to-continuum ratio. The latter is calculated dividing 
the velocity integrated line-to-continuum ratios by the linewidth. The errors given are formal standard deviations that do not 
include calibration uncertainties. Col. 8: Continuum flux density of image A w.r.t. the total continuum flux density. The error is 
about ±0.03. Col. 9: Peak "apparent" optical depth, i.e. r = -ln(l-|TL|/T c ), assuming uniform coverage of image A and an H2CO 
excitation temperature that is small w.r.t. the continuum background. 



systematically decreasing flux densities beyond ^15 GHz, 
our ^85 GHz fluxes are consistent with the correspond- 
ing part of the continuum spectrum shown by Combes & 
Wiklind lITWTj) and Combes et al. (jTWF|) that suggest a 
spectral index of a ~ -0.25 (S v oc v a ). The minimum in 
our "spectrum" at 43.2 GHz (Tabled]) may be caused by 
time variability or by a substantial amount of missing flux. 
Among the data presented here, the 43.2 GHz Q-band map 
(Fig. 2) has the highest angular resolution so that all spatial 
components except images A and B are resolved out. 

This might also affect our measured A/B flux density 
ratios. While B is located near the center of the Einstein 
ring, A is placed near its northeastern edge. Thus estimated 
flux densities for B may be contaminated by Einstein ring 
emission at 8.6 and 14.1 GHz, but not at 43.2 GHz. This 
would be consistent with lower ratios at the two lower fre- 
quencies. Nevertheless, our A/B flux density ratios, ~ 3.3 
(at 8.6 and 14.1GHz) and ~3.5 (at 43.2 GHz) are, within 
the formal errors (see Table [2]), consistent with a single 
value. Furthermore, they also lie within the range of values 
obtained by Mittal et al. (pZOOOl their Table 3 and Fig. 5b) 
for frequencies 8.4-15.4 GHz. Thus we confirm that there 
is no evidence for a significant amount of free-free absorp- 
tion between 8.6 and 43.2 GHz (see Mittal et al. 120071 for a 
model of such an ionized absorbing layer). 

4,2. Spectral analysis 
4.2.1. General aspects 

Table [3] gives integrated (Col. 4) and peak (Col. 7) line-to- 
continuum ratios of the six lines shown in Fig. [3] as well as 
the fractions of the total continuum flux originating from 
image A (Col. 8) and the resulting peak apparent optical 
depths (Col. 9). Assuming that only image A is absorbed 
by a foreground molecular cloud (see Sects. 1 and 3.1 and 
Mittal et al. 2007) and that the coverage is uniform, the 
rotational H2CO transitions have optical depths 0.3 < t < 
0.9, while the K-doublet opacities are at least an order of 
magnitude lower. 



Henkel et al. (2005) interpreted the trend of increas- 
ing optical depth of molecular lines with rising frequency 
in terms of a uniform molecular absorber located toward 
the line-of-sight to the peak of a radio continuum source. 
In this scenario, the size of the continuum source decreases 
with increasing frequency. In view of the trend in apparent 
optical depth that is obvious when comparing lines at dm-, 
cm- and mm-wavelengths (e.g., Carilli et al. 119931 Menten 
& ReidH996l Wiklind & Combes [[995]) it may seem to be 
surprising that the J=2-l para-H 2 CO line, observed at a 
similar frequency as the corresponding ortho-lines is less 
optically thick (Table [3]). The difference may be caused by 
a high ortho- to para-H2CO abundance ratio or by a low 
excitation of the absorbing molecular gas (see Combes et 
al. fTMTI for CS exc itation and Combes & Wiklind fl997l and 
Gerin et al. 119971 for limits to the H 2 and 13 CO excita- 
tion). For ortho-H2CO, the J=2-l lines are the lowest rota- 
tional transitions. Para-H 2 CO, however, also has a J=l-0 
line. Thus the J=2-l para-H 2 CO line is not the lowest ro- 
tational line in its respective K a ladder (Fig.[l| and Table[3] 
shows indeed that the J=l-0 line has an optical depth that 
is similar to that of the J=2-l line, in spite of the lower sta- 
tistical weights of its two rotational states. However, it has 
to be emphasized that the assumption of a uniform molec- 
ular cloud is highly unrealistic. The use of this assumption 
is only justified by the limited spatial resolution available 
when studying the molecular content of a galaxy at redshift 
z=0.68. 

4.2.2. Continuum source morphology 

Sensitive high resolution continuum maps of image A at 
8.4 GHz show a source of size ~10xl0mas 2 that is edge- 
brightened on its south-western side (Fig. 2 of Biggs et 
al. I2003|) . Tangential stretching by a factor of 3-4 along 
position angle P.A. ~ -30° almost matches the elonga- 
tion caused by a jet at P.A. ~ +60°. At low frequencies 
(<8GHz), the morphology of component A may be af- 
fected by free-free absorption (Mittal et al. l2007[) . At higher 
frequencies, VLBI continuum measurements of B0218+357 
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have been made up to 22 and 43 GHz (Porcas & Patnaik 
1996; Porcas 2004), but sensitivities are not high enough 
to reveal additional spatial fine structure. It would not be 
a surprise if the jet contribution would fade at higher fre- 
quencies leaving the tangentially stretched core that may 
coincide in position with a spiral arm of the lensing galaxy 
(see Fig. 3 of York et al. 2005). Assuming that the core 
shows a brightness temperature and flux density not dras- 
tically varying with frequency in the 10-100 GHz frequency 
range (e.g., Blandford & Konigl iT979p . its solid angle would 
vary as SI oc j/" 2 ---" 1 (see also Lobanov 1998) between cm 
and mm- wavelengths. With a spectral index of a ~ -0.25 
(Sect. 4.1) the mm-wave flux of B0218+357 is about half of 
that at cm- wavelengths, which is compatible with a fading 
jet and a flat spectrum core. 

4.2.3. LVG simulations 

To test the qualitative scenario outlined in Sect. 4.2.1 
and to account for the background continuum morphol- 
ogy sketched in Sect. 4.2.2, we use Large Velocity Gradient 
(LVG) models for para- and ortho-H 2 CO, accounting for 
the 41 and 40 lowest energy levels of each species up to 
300 cm" 1 above the ground state. We adopt collision rates 
from Green (|1991|) and assume a spherical cloud geometry. 
The choice of a particular cloud geometry can affect the re- 
sulting densities up to half an order of magnitude, but only 
if the lines are optically thick. Applying a plane-parallel in- 
stead of a spherical cloud geometry, resulting particle den- 
sities could be lower by up to this amount. 

Since line shapes are identical within the limits of noise, 
we can try to reproduce the various peak optical depths 
of the main velocity component displayed in Table [3] As 
we shell see, main constraints are obtained with the four 
ortho-H2CO transitions, while the para-H2CO lines serve 
as a measure of ortho-to-para abundance ratios. 

The simplest case, ignoring collisions: 

As a first approach, we neglect collisional excitation as well 
as frequency dependent variations of the underlying contin- 
uum morphology (Sect. 4.2.2). Since the data do not signif- 
icantly constrain the temperature of the cosmic microwave 
background (CMB), we adopt T CMB = 2.73 x (l+z)K = 
4.60 K = T cx (T ex : excitation temperature). Table 2] displays 
the resulting optical depths (Model I) with minimized re- 
duced x 2 values of 7.1 for the four ortho- and 1.7 for the two 
para-H2CO transitions (x 2 /(N — P), N: number of lines, P: 
number of free parameters; P rtho = 2 (density and column 
density), P para = 1 (column density)). The total column 
density is V(H 2 CO) ~ 3.9 xlO 13 cm" 2 and the ortho to para 
abundance ratio becomes ~2.2. Modeled optical depths of 
the two para-H2CO transitions are well within the 2a limits 
of the observational errors. For the two cm- wave K-doublet 
ortho-lines (see Fig. [I]), however, modeled opacities appear 
to be too high, while for the mm-wave J=2— 1 lines modeled 
optical depths tend to be too low. Therefore the calculated 
ortho- values raise some doubts about the quality of the fit. 

Introducing collisional processes: 

Toward B0218+357, neglecting collisions is a potentially 
realistic assumption, when analyzing diatomic or linear 
molecules with a CO-like energy level scheme (see Combes 
& WiklindHMll Gerin et al. fLWl and Henkel et al. [2005)1 . 



5 

Table 4. Optical depths a) 



Line 


U bserved 


Model 


Model 






I 




i,„_i,, (0) 


09*1+0 007 

U. KJZiO -l—VJ .UU 1 


0.028 


0.036 


2n — 2i2 (O) 


0.006±0.002 


0.010 


0.007 


212-I11 (O) 


0.493±0.090 


0.436 


0.183 


2u-lio (O) 


0.890±0.151 


0.426 


0.160 


loi-Ooo (P) 


0.331±0.068 


0.257 


0.250 


202-I01 (P) 


0.311±0.057 


0.352 


0.346 


'2 

Xreduced, ortho 




7.1 


18.9 


Xrcduccd,para 




1.7 


1.8 


V(H 2 CO) 




3.9 


2.6 


Ortho/Para 




2.2 


0.9 



a) Col. 1: H2CO transition with O and P denoting ortho- and 
para-H2CO, respectively. The upper two lines refer to K-doublet 
transitions. In the lower part of the Table, reduced x 2 values, 
total H2CO column densities from the +7kms _1 component in 
units of 10 13 cm -2 , and H2CO ortho-to-para abundance ratios 
are given. 

Col. 2: Observed apparent optical depths taken from the last 
column of Table [3] Given errors are larger than those in Table [3j 
also accounting for a calibration uncertainty of ±15%. 
Col. 3: LVG results exclusively accounting for microwave back- 
ground excitation (Tcmb = 4.6 K = T cx ; Model I). 
Col. 4: Same as Col. 3, but including collisional excitation 
(n(H 2 ) = 500 cm" 3 ; Model II). For details, see Sect. 4.2.3. 



H 2 CO, however, critically responds to collisional excitation 
even at very low densities. This is a consequence of the col- 
lisional cooling of the K-doublet lines (see e.g., Garrison et 
al. U975p that sets in at densities as low as n(H2)^100 cm" 3 , 
while collisional heating of the rotational transitions is still 
far from being significant. Collisional cooling of the lio — In 
and 2n — 2i2 lines is rather independent of the chosen ki- 
netic temperature as long as T^ n > 10 K. Here we adopt 
Tkin = 55 K (Henkel et al. I2005[) . Reduced T ex values cause 
increased optical depths in the K-doublet lines. As a conse- 
quence, the optical depths of the lio — In K-doublct and the 
rotational ortho-H2CO lines become even more similar than 
in model I and a reasonable fit cannot be obtained. This is 
exemplified for a density of n(H2) = 500 cm" 3 in Model II of 
Tabled] which is characterized by a reduced x 2 value of 18.9 
for the ortho-transitions. Higher densities, i.e. 500 cm" 3 < 
rt(H2) < 10 5 cm" 3 , pose similar problems, while densities 
JJlO 5 cm" 3 severely contradict additional constraints that 
have been discussed by Combes et al. (1997), Gerin et al. 
(fT997j) . and Henkel et al. (|2005|) . 

A frequency dependent radio continuum morphology: 
In view of the problems outlined above, the best way to 
reconcile observational and model data is to account for 
the frequency dependence of the continuum morphology 
(Sect. 4.2.2). A background continuum source covering fac- 
tor / c , increasing with frequency, would (1) reduce the dif- 
ferences in optical depths between the lio — 1 11 K-doublet 
and J=2-l rotational lines and would (2) raise the optical 
depth ratio of the K-doublet lines. At 137 GHz, / c ~l for 
image A (Wiklind & Combes HHMD- At 0.85 GHz, we get 
from the redshifted 1.7 GHz OH absorption lines and the 
continuum flux of the entire source (Kanekar et al. I2003|) a 
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Table 5. Optical depths accounting for source coverage 31 ' 



Line 


Observed 


Model 


Model 


Model 






III 


IV 


v 


±10 ^11 V / 


n 1 1 3+0 039 


0.046 


0.086 


0.120 


2n — 2i9 (0) 


0.016±0.006 


0.017 


0.020 


0.024 


212-I11 (0) 


0.539±0.098 


0.726 


0.641 


0.610 


2u-lio (0) 


0.941±0.160 


0.710 


0.591 


0.533 


I01-O00 (P) 


0.503±0.103 


0.286 


0.285 


0.283 


202-I01 (P) 


0.335±0.061 


0.393 


0.393 


0.392 


Xreduced, ortho 




4.9 


3.5 


4.4 


Xr educed, par a 




5.3 


5.4 


5.4 


iV(H 2 CO) 




5.4 


5.2 


4.8 


Ortho/Para 




3.0 


2.8 


2.5 



a) Col. 1: H 2 CO transition with O and P denoting ortho- and 
para-H 2 CO, respectively. The upper two lines refer to K-doublet 
transitions. In the lower part of the Table, reduced x 2 values, 
total H2CO column densities from the +7kms~ 1 component in 
units of 10 13 cm -2 , and H2CO ortho-to-para abundance ratios 
are given. 

Col. 2: Observed apparent optical depths (see Table Q, but ac- 
counting for a frequency dependence of the continuum source 
covering factor, f c — (z^ghz/100) ' 5 . 

Col. 3: LVG results exclusively accounting for microwave back- 
ground excitation (T C mb = 4.6 K = T ex ; Model III). 
Cols. 4 and 5: Same as Col. 3, but including collisional excitation 
with n(H 2 ) = 200 and 500 cm -3 , respectively (Models IV and 
V). For details, see Sect. 4.2.3. 



peak apparent optical depth of 0.01, implying that the cov- 
ering factor f c is <;0.01. Assuming that f c — 1 is reached at 
exactly 100 GHz, this corresponds to f c oc j/ - 97 (choosing 
instead 70 or 130 GHz would yield similar exponents, i.e. 
1.04 or 0.91). The lower f c limit at 0.85 GHz is, however, too 
small. Firstly, the 1.667 and 1.665 GHz OH lines show the 
intensity ratio expected in the optically thin limit at local 
thermodynamic equilibrium conditions (1.8:1). Therefore, 
unlike for CO, the OH optical depths are well below unity 
(t.^0.5). This increases the attainable line-to-continuum ra- 
tio in the optically thick case (i.e. the source covering fac- 
tor) at 0.85 GHz by at least a factor of two and yields f c oc 
^0.0. ..0.8 Secondly, we also have to account for the fact that 
at low frequencies the Einstein ring becomes prominent and 
that image A is not dominating the total flux density (for 
free-free absorption of image A, see Mittal et al. l2006|) . Only 
taking the continuum flux from image A thus further raises 
the 0.85 GHz line to continuum ratio and f c cx v a - a ■ °- 5 be- 
comes a more realistic estimate. This is a weaker frequency 
dependence than the hypothesized one for fi, the contin- 
uum solid angle (Sect. 4.2.2), but agrees with the observed 
overall spectral index of the continuum radiation. 

In order to study the case with an exponent for f c of 0.5 
that differs most from a scenario with a frequency indepen- 
dent f c value, Table[H]displays the / c -corrected observed op- 
tical depths as well as LVG model results for densities n(H2) 
= 0, 200, and 500cm" 3 . Excluding collisions (Model III), 
the two K-doublet lines show optical depths that are too 
similar for a good fit. Minimized reduced % 2 values are 4.9 
and 5.3 for the ortho and para-H2CO lines, respectively. 
A density of 200 cm -3 (Model IV) reproduces much bet- 
ter the K-doublet lines and shows with reduced \ 2 val- 
ues of 3.5 and 5.4 the highest probablity. Total H2CO col- 



umn densities and ortho-to-para abundance ratios become 
5.2xl0 13 cm~ 2 and 2.8. At higher densities we get the same 
effect as in the case of a frequency independent f c value. 
Opacities of K-doublet and rotational ortho-H2CO transi- 
tions slowly converge and \ 2 values gradually increase with 
rising density. Model V shows this for n(Ha) = 500 cm~ 3 , 
while for n(H2) = 1000 and 3000cm -3 , reduced x 2 values 
for the ortho-H2CO lines become 6.0 and 8.8, respectively. 
The significant difference in the measured optical depths of 
the two J=2— 1 lines of ortho-H2CO is not reproduced in 
any model and deserves observational confirmation. 

From Table [5] it is apparent that the fits to the para- 
H2CO lines are almost independent of density, because ro- 
tational excitation is minimal at the densities considered 
here. The LVG simulations of the para-H2CO lines in mod- 
els III, IV, and V are not as good as in models I and II. 
While obviously not being suitable for a reliable Tomb de- 
termination, calculated opacities arc still within the 1 and 
3cr limits of the observed values for the 2o2~ loi & n d loi - 0oo 
lines, respectively. 

4.2.4. More general implications 

To summarize, the successful models (I, III, V, and par- 
ticularly IV) agree well with each other, both with re- 
spect to column density and ortho-to-para abundance ratio. 
The density of the gas is n(H2) < 10 3 cm~ 3 , most likely 
~200cm -3 , supporting the view that the cloud is diffuse 
and has peculiar physical parameters rarely found in the 
nearby universe (cf., Henkel et al. I2005p . For the best sce- 
nario, Tki n = 55 K (from NH 3 ) and n(H 2 ) = 200 cm~ 3 (from 
H 2 CO), we have calculated excitation temperatures, assum- 
ing a cosmic background radiation of 4.6 K. In the optically 
thin case, excitation temperatures of the three ground rota- 
tional transitions of CS, HCN, HNC, HCO+, and N 2 H+ are 
within 1% of the assumed CMB background temperature 
and should thus be extremely useful for future measure- 
ments aiming at 2cmb as a function of redshift. Rare CO 
isotopomers, however, are not as useful. 

There is a second velocity component, at -7kms _1 (see 
Fig. [3] and Table [3]), that is also apparent (mostly in the 
form of skewed profiles) in transitions of other molecules 
(Wiklind fc Com bes ITM51 Com bes et al. [TM71 Combes 
& Wiklind IT9981 Kanekar et al. [2003|) . The broad com- 
ponent seen in ammonia (NH3) at a heliocentric velocity 
of ~0kms _1 (Henkel et al. 120051 use Local Standard of 
Rest (LSR) velocities, while here heliocentric velocities are 
given) is a superposition of hyperfine satellites of the ~ 
+6kms~ 1 and the entire -7kms _1 line and may thus serve 
as a rough measure of the relative importance of the two 
kinematic NH3 features. As is the case for ammonia, the 
-7kms _1 component of formaldehyde is too weak for de- 
tailed modeling but adds ~15% to the resulting column 
densities given in Tables 2] and 

The ortho-H 2 CO abundance appears to be slightly 
larger than that estimated by Menten & Reid (1996). 
Having carried out multiline studies of ammonia (NH3) 
and formaldehyde toward B0218+357, we can now compare 
their abundances. For the total column densities (because 
of hyperfine splitting, the components are not well sepa- 
rated in NH 3 ) we fi nd N(K 2 CO) - 0.6 x N(NH 3 ) (for NH 3 , 
Henkel et al. 12005]) . The A(H 2 CO)/7V(NH 3 ) abundance ra- 
tio is four times smaller than the average ratio estimated 
for local diffuse clouds observed in absorption (Liszt et al. 
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2006). In view of a possible overabundance of the more 
'primary' 12 C and 16 nuclei relative to the mainly 'sec- 
ondary' 14 N nucleus in spiral galaxies of the distant past 
(see e.g., Wheeler et al. 1989 for a review on metallicity 
dependent CNO abundances) this may seem to be surpris- 
ing. However, the galactic H 2 CO abundances are based on 
two H2CO transitions (versus six here), while NH3 data are 
from measurements of two spectral lines only (versus four in 
the case of B0218+357). For only a few diffuse clouds galac- 
tic NH3 absorption profiles have been measured. Therefore 
the 'local' average ratio (see also Nash I1990|) may be less 
certain than the ratio presented here for a single distant 
gravitational lens alone. While our H2CO column density 
is 2.5 times higher than the highest column density among 
the sources summarized by Liszt et al. (|2006p . our absorp- 
tion profile is wider by at least this amount, thus yielding 
a similar column density per kms -1 velocity interval. If 
H2CO is preferentially formed on dust grains while NH3 is 
arising in the gas phase from either N + or N (e.g., Liszt 
et al. I2006P , a good correlation between the abundances of 
the two molecules may not be expected. 

An ortho-to-para abundance ratio of almost three and a 
kinetic temperature of 55 K (Henkel et al. I2005[) are consis- 
tent with the correlation between kinetic temperature and 
ortho-to-para abundance ratio shown in Fig. 10 of Kahanc 
et al. (|1984p . Whether this indicates gas phase H2CO for- 
mation in a warm medium or whether adsorption and des- 
orption on grain surfaces can lead to efficient thermalization 
(Mangum & Wootten ll993j) remains to be seen. 



5. Conclusions 

The detection of six H2CO absorption lines toward the 
gravitational lens B0218+357 reveals surprisingly detailed 
information on molecular gas at redshift z=0.68. Our anal- 
ysis leads to the following main results: 

— The continuum flux density ratio of component A to 
B is ~3.3 between 8.6 and 43 GHz, does not vary sig- 
nificantly, and is thus not severely affected by selective 
free-free absorption in this frequency range. 

— The loi — Ooo line of para-H2CO is mainly absorbing the 
radio continuum from image A, consistent with previ- 
ously published results from the 2n — 2i 2 line of ortho- 
H 2 CO. 

— The H 2 CO column density of the main spectral fea- 
ture is ~5xl0 13 cm -2 . The mm-wave lines show a sec- 
ond, blue-shifted velocity component that is displaced 
by about lSkms^ 1 from the main component. It con- 
tributes another ^15% to the total column density. 

— The best fits to the ortho-H 2 CO molecular lines involve 
the assumption of a continuum cloud coverage that, be- 
tween 0.85 and 100 GHz, is likely frequency dependent 
•:/, oc iAQ-0-5). 

— The gas density is n(H2) < 1000cm 3 , most likely 
^200 cm -3 , confirming the view that the absorption is 
arising from a diffuse molecular environment. 

— With Tkin and n(H2) known, optically thin lines of the 
three ground rotational transitions of CS, HCN, HNC, 
HCO + , and N2H 4 " (not CO) will be ideal measures of 
the cosmic microwave background at z=0.68. 

— The total H 2 CO column density is about 60% of that 
obtained for NH3. This is four times below the aver- 



age ratio deduced from local galactic diffuse gas seen in 
absorption. 

— The ortho-to-para H2CO abundance ratio is large (2.0- 
3.0), reflecting the kinetic temperature of the cloud. 
Whether this is caused by cloud formation in a warm 
medium or whether this is due to efficient thermaliza- 
tion on dust grain mantles remains open. 
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Fig. 3. H 2 CO lines observed toward B0218+357 and 
Gaussian fits. Given are line-to-continuum ratios (the en- 
tire continuum flux of B0218+357 was taken) as a function 
of heliocentric velocity adopting a redshift of z—0. 68466. 
The labels in the lower left corner of each spectrum show 
transition and rest frequency. O and P refer to ortho- and 
para-H2CO, respectively. 



x — offset (arcsec) 



Fig. 2. VLA A-array continuum images of B0218+357 
taken at Q-, U-, and X-band (top to bottom). In each im- 
age, the lowest contour line represents four times the rms 
noise level, which is a — 0.21, 2.1, and 0.5 mJy beam -1 
for the Q-, U, and X-band images, respectively. Contour 
levels increase by a factor of 4 for the Q-band image and 
by a factor of 2 for the U- and the X-band images. Peak 
flux densities are 0.267, 0.99, and 1.0 Jy beam -1 . The full 
width to half maximum dimensions of the restoring beam 
are shown in the lower left corner of each panel. The origin 
of each image is the phase center position given in Sect. 2.1. 



